AD 697 697 


A STATISTICAL TREATMENT OF VARIOUS CLASSES OF 
GUNNERY ERRORS AND THE CALCULATION OF HIT 
PROBABILITY 


Hermann J. Helgert 


Cornell Aeronautical Laboratory, Incorporated 
Buffalo, New York 


30 September 1969 


Distributed ... ‘to foster, serve 
and promote the nation's 
economic development 

and technological 

advancement.’ 


CLEARINGHOUSE 


FOR FEDERAL SCIENTIFIC AND TECHNICAL INFORMATION 


U.S. DEPARTMENT OF COMMERCE/ National Bureau of Standards 


This document has been approved for public release and $ate.~ . Ave fe 


¥e.c21e Copy 


fen SARS YORE IT OR ae PEPER PES AE AT ne OV OL BOSS SO EL SERIE PEI TS piers mee mens 


CLEARED 
FROPEN PUBLICATION 


19 NOV 1969 


Chief of Naval Material (PA} 


A STATISTICAL TREATMENT OF VARIOUS CLASSES OF 
GUNNERY ERRORS AND THE CALCULATION OF HIT PROBABILITY 


AD 697697 


PROJECT SECANT Il 
FINAL REPORT 
VOLUME I 


30 September 1969 


st 7 


CAL Report No. GM-2717-D-2 


Hermann Josef Helgert (| 2 : 


DEC 5 = 1969 


Lyon UW LD 


Cc 


PREPARED UNDER CONTRACT N00178-69-C-0034 FOR THE 
U.S. NAVAL WEAPONS LABORATORY, DAHLGREN, VIRGINIA 
BY CORNELL AERONAUTICAL LABORATORY, INC., BUFFALO, NEW YORK 
Reproduced by the 
CLEARINGHOUSE 
for Fedwral Sciantitic & Technical 
Information Springlield Va. 22151 


DISTRIBUTION OF THIS NCCUMIENT IS CALENITED, 


t> 
co 


Best Available Copy 


Smee. 


TCGESSION OO 
WHE Sos, 


cesT 

pot purr Sees Ty 
BNANNOUNCED cmt 
JUSTIFICATION eer > 


_ 
saa | 


BY Sccagetnn 
DISTRIBUTION AVAILABILITY GuNEe 


~pist. 


L gyANL. ane & SPEC aL 


~*~ 
— may 


A STATISTICAL TREATMENT OF VARIOUS CLASSES OF 
GUNNERY ERRORS AND THE CALCULATION OF HIT PROBABILITY 


PROJECT SECANT II 
FINAL REPORT 
VOLUME I 


30 September 1969 


CAL Report No. GM-2717-D-2 


by: 
Hermann Josef Helgert 


PREPARED UNDER CONTRACT N00178-69-C-0034 FOR THE 
U.S. NAVAL WEAPONS LABORATORY, DAHLGREN, VIRGINIA 
BY CORNELL AERONAUTICAL LABORATORY, INC., BUFFALO, NEW YORK 


NOTICE TO USERS 


Portions of this document have been judged by the Clearinghouse 
to be of poor reproduction quality and not fully legible. However, inan 
effort to make as much information as possible available to the public, 
the Clearinghouse sells this document with the understanding that if the 
user is not satisfied, the document may be returned for refund. 


If you return this document, please include this rotice together with 
the IBM order card (label) to: 


Clearinghouse 
Attn: 152.12 
Springfield, Va. 2215] 


‘ 
aa TABLE OF CONTENTS 
Section Title Page Ne. 
ACKNOWLEDGEMENTS G4. i web dae oe hee Ae es eta ae AY 
ABSTRACT . bs fo Vv 
1 INTRODUCTION AND SUMMARY... 2... eee ee ee te ee 1 
2 MISS DISTANCE STATISTICS . 3 
3 COMPUTING THE MISS DISTANCE STATISTICS FROM FIRING DATA . 11 
4 DELIVERY: ACCURACY: uclie 2 ee ait Bae ae Ste a! OG 
5 SINGLE SHOT HIT PROBABILITY, CUMULATIVE HIT PROBABILITY, AND 
EXPECTED NUMBER: OF HITS: “a. 2.97 fo uis Seo ek ae eee ORS 
6 APDETATUED: EXAMPLE % cy G04: 5a daw Sse Se. gh Be ge alg. 28 
; 7 TABLES. “AND -<GRAPHS. ue: Se ac We Ra ee oe ee a Ve “80 
8 REPERENGES 3. e- ai si @ Sg ne Oe WE ee a ae re, SS 
APPENDIX ‘TABLE OF CUMULATIVE HIT PROBABILITY .......... 4 Ack 
poe iii 


ACKNOWLEDGEMENTS 


The author gratefully acknowledges the cooperation and assistance 
provided by Dr. W. A. Kemper of the U.S. Naval Weapons Laboratory, Dahlgren, 
Virginia. In addition, suggestions and assistance were provided by Fred J. Tanis 
and John J. Winkleman, programming for the IBM 360/65 computer by Steven L. 


Sackett, and final report review and editing by Fdward J.H. Lane, all of the CAL 
staff. 


ABSTRACT 


The work reported in this final report was performed by Cornell 
Aeronautical Laboratory, Inc. for the U. S. Naval Weapons Laboratory, Dahlgren, 
Virginia, under contract No. N00178-69-L-0034 and is submitted in fulfillment 


of that contract. 


Volume I of the report is a manual and presents, in simplified form, 


methods and exampies of computing the hit capability of Naval Gun Systems 


{ 
i 


through the use of tables, graphs and equations. The methods permit considera- 
tion of any correlations that may exist between successive rounds in a salvo. 
This voiume, however, does not treat the problem of establishing the degree of 


correlation. 


Although the equations derived are general in nature, the tabular data 
appended provide estimates of achieving at least one hit per salvo in the 
special situation where the range and cross-range components of the error 
statistics are eyual and the turgets are nearly square or circular. There are 
techniques whereby non-cqual components of error statistics may be converted 
to equivalent equal components; however,their applicability to the problem at 
hand has not been examined and is considered beyond the scope of this study. 


Single shot hit probabilities may be derived from the case where the salvo 
size equals one. Expected number of hits on the target per salvo may be 
acquired by multiplying the single shot hit probability by the number of shots 
in the salvo,assuming there is no reason to belicve that the single shot hit 
probability changes during the salvo. Where this is not the case the salvo may 
be subdivided into intervals in which the single shot probability is essentially 
constant, the expected number of hits computed for each interval and ths values 
totaled. 


‘ Volume II of the report defines a program of analysis and testing which would 
provide, for any situation, estimates of the component error parameters required 
! to enter the tabular data of Volume |. 
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1. INTRODUCTION AND SUMMARY 


The ultimate worth of a gun system lies in its ability to destroy the target. 
The effectiveness with which it accomplishes this task is thus a fundamental 
parameter of the system. For a number of reasons it is important to have good 
estimates of weapon effectiveness; they can, for example, provide a basis for 
competitive assessment of the gun system against other existing or planned 
weapon systems and are a requirement for proper logistic and firepower mission 
planning. 


Weapon effectiveness is influenced by three quantities: the target char- 
acteristics or vulnerability, the type of warhead and lethal mechanism used, and 


the accuracy with which the weapon can deliver the warhead on or near the target. 


For hard targets the interaction of the first two factors is defined in 
such terms as target vulnerable area or the probability of a kill per hit. For 
soft targets which may he destroyed or damaged by near misses of blast-frag- 
mentation type warheads, this interaction is usually defined in terms of warhead 
lethal area. The vulnerability of both target types then, is specified in terms 
which recognize the type of warhead under consideration, and the latter is 
itself determined by the affects to be achieved at the target. The problem we 
are concerned with in this report, however, is the computation of the third 
factor, delivery accuracy. 


While the solution to this problem, under certain assumptions, is known and 
documented, reference 1, the numerical application of these results is fraught 
with difficulties. First the solution draws heavily on the theory of discrete, 


stcchastic processes and requires for its understanding at least some acquain- 
tance with this discipline. Second, the mathematical format of the solution 


is complicated and does not lend itself easily to numerical computation. 


For these reasons it is desirable to present the problem and its solution 
in a somewhat simplified form in order to make the results available to a wider 
class of users who would otherwise find them intractable. A previous report, 
reference 2, treated this problem to a large degree, but lacked so.ewhat in 
applicability. The objective of this manual is to present a more comprehen- 


sive treatment which avoids this limitation. 


The computation of delivery accuracy requires, first of all, knowledre of 
the statistical pronerties of the proiectile miss distances. The nature of 
these quantities and the basic assumptions inherent in the solution are discussed 
in Section 2. 


Section 3 treats the problem of estimating the miss distance statistics 


from records of actua! test firings. 


Messures of delivery accuracy and the computational formulas are presented 
in the following section. Three of the most important measures are discussed 
separately in Section S. The manual concludes with a detailed example of the 
computations involved, and a graphical and numerical tabulation of cumulative 


hit probability. 
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2, MISS DISTANCE STATISTICS 


Consider a pun firing projectiles at a surface target. Because of 
certain errors inherent in the gun system as a whole the observed impact 
point will, in g»neral, differ from the intended one by an amount which is 
randon from shot to shot. This difference is a composite of three funda- 


mental errors (see Figure 1): 


1. Systematic errors--those which can be considered constant for the 
duration of an engagement but which vary from one engagement to 


the next. 


2. Time-varying ecrrors--those which vary significantly during the 
engagement but whose rate of variation is slow compared to the 


firing rate. 


3. kound-to-round dispersion errors--those which vary in an uncorre- 


lated manner from one round to the next. 


Systematic crrors arise from many sources such as imbalances in the 
servos and wrong pain scttings in the amplifiers of the radar and fire con- 
trol computer, inaccuracies in wind, temperature and air density estimates, 
failure to fully account for changes in initial round velocity due to barrel 
wear, navigational errors, crrors in locating the target in indirect fire, 
etc. Their net effect is to impart a bias on the center of impact points, 
whose value is constant but unknown for a particular cngapement and varics 


from one cnpagement to the next. 


The time varying crrors are duc to such thinrs as changes in tube sag 
because of heating effects, ship flexure, incorrect radar tracking data which 
is processed through the fire control computer, and chanpes in metcorological 
factors. These give rise to aim-wander, a term which derives from the fact 
that the path traced by the intersection of the pun mean line-of-sight and 
a plane perpendicular to it would, as a function of time, appear to be wan- 
dering in a more or less random fashion (see Figure 2). As the gun fires a 


succession of rounds, it samples the aim-wander path at the instants of 


firing. 
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Figure 2 TYPICAL AIM-WANDER PATH 


Finally, round-to-round dispersion errors are caused by the combined effect 

of round-to-round variation in shell manufacture, powder weight, moisture 

: content and temperature, and short-term turbulence in the state of the atmos- 
phere during the time of firing. In this list must be included a factor peculiar 
to guns mounted on unstable platforms such as ships, namely round-to-round 
variations in range and deflection caused by translational motion of the gun 
barrel at the instant of firing, and th: aature of the recoil of guns mounted 
on such unstable platforms. 


Table 1 contains a listing of the important error sources in a typical 
weapon system such as the 5"/54 cal. gun with Mk 68 fire control system together 
with a classification of the impact point errors they give rise to. Not all of 
these sources, however, are fi.ctors in all modes of fire. Note that this 
classification is not unique in the sense that some of t:.. system errors may 
cause more than one type of impact error. Neither are the classifications 
indicated in the table to be treated as final. At this time the teble should 
be looked upon as being illustrative of the fact that the error is considered 
in this tentative assessment to be predominantly of that type. The classifica- 
tion or these error sources will receive more comprehensive treatment in Vol. II. 
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TABLE 1. ERROR CLASSIFICATION 


syst. [TV] Disp. 


Gun Director 
Target Range 
Target Bearing 
Target Altitude 
Director Misalignment in Direct Fire 
Radar Tracking Errors and Lags 


Navigation System 
Ship Position 
Errors in Locating Target Position and Vel. 
Spotter Errors 


Gyro Compass 
Ship Heading 


Pitometer Log 
Ship Speed 


Offset Inputs 
Range 
Deflection 


Stable Element 
Stabilization Errors 


Meteorological Data 
Wind Speed and Direction 
Air Temperature and Density 


TABLE 1. (CONT) 


[ Syst. | T-v | Disp. | 


8. Fire Control Computer 
Target Prediction Errors 
Transients 
Blunders in Operating the Comriter 
Amplifier Imbalances and Lags 


Coordinate Conversions 


9. Data Transmission 


Synchro Errors 


10. Alignment 
Roll Path 
Parallax 
Ship Flexure and Bending 


11. Weapon Factors 
Tube Sag 
Tube Heating 
Jump Variation 
Initial Velocity Variations 
Drag Variations 
Gun Wear 


12. Weapon Response 
Response to Gun Orders 
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It has, in the past, been assumed that the time-varying error is either 
absent or of the same nature as the round-to-round dispension error. There 
may be cases, however, when this assumption does not hold true. Computation 
of weapon effectiveness then may be grossly in error. The significance of the 
time-varying error lies in the fact that its values at different firing 
instances may be correlated. [It is only in the case where this correlation 
is negligibly small that the above simplification may safely be made. 


All subsequent work concerning the evaluation of gun fire effective- 
ness is based on two assumptions: 


(a) The three types of errors are Gaussianly distributed and independent 
of each other. 


(b) That statistical nature of the errors does not change with time. 


The assumption of a Gaussian distribution has been well established 
for the round-to-round dispersion errors, but is somewhat questionable for 
the other types of errors. 


Assumption (b) requires, among other things, a stationary gun-target 
geometry, i.e., no change in range and line of fire, no spotting adjustment 
of the aim point during the firing and constant meterological conditions. 
Usually the engagements are short enough so that changes in the meteorological 
conditions may be discounted. 


If there is appreciable motion between gun and target during the 
firing, the engagement can be split into sections during each of which the 
relative motion is negligible. Another device which is sometimes helpful is 
to define a meai. coordinate system in the sense that the excursions of the 
line of fire and range from this system are small during the firing. 


If,during the course of firing a salvo being analyzed by technique: of 
Sec. 6, spotting corrections are introduced into the aim, the magnitude of 
these corrections must be recorded. They can then be subtracted from the impact 
points and the resulting data will be statistically stationary. 


If (a) and (b) are satisfied, the three types of errors are uniquely 
determined by the following statistical parameters, measured in an x-y coor- 
dinate system located on the target (see Figure 1): 


Someone se : ——— 


Systematic Error - standard deviations s. and Sy over all 
possible engagements and actual values 
u, and uy, on one engagement 

Time-varying Error - standard deviations a, and a) and 
correlation functions 2. (j) and 


Ay, (i) 


Round-to-Round - Standard deviations d. and J, 


Dispersion krrors 


The standard deviations are the major and minor diameters of an ellipse cen- 
tered at the origin of the error vector of Figure 1, within which approxi- 
mately 40% of the errors fall. 


The correlation functions are a measure of the statistical con- 
straint placed by the time-varying errors on succeeding values of these 
errors, For a constant firing rate, — depends only on the difference in 


round nunber and typically takes the form indicated in Figure 3. 


The quantity j is the time difference between particular rounds 
and A(j) is the correlation in the time-varying errors for these rounds. 
Thus, for example Y (1) is the correlation in these errors for any two rounds 
one unit of time apart, (2) the same quantity for rounds two units apart, 
etc. Note that 70) is the correlation between an error and itself. 


In many applications, the functions A are too complicated to be of 
general use. In such cases, these functions are replaced by constants rT. and 
r_, chosen so that the delivery accuracy computed on the basis of these con- 
scant correlations is a close approximation to the true value. Reference (1) 
indicates that in situations where the correlation constants exceed approxi- 
mately 0.2, the hit probabilities would be significantly different from those 
assuming round to round independence. The size of this correlation constant 
depends on the nature of the correlation function. If may be possible to 
characterize this function once the various contributing processes to the 
time varying error have been thoroughly examined in a program of the type 
discussed in Volume II. The physical interpretation of this change is that 
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we replace a system in which the shots are fired serially by an equivalent 


system which fires its projectiles in « simultaneous burst. 


The parameters uy and u_ define the true average aim point on a fixed 


engagement and represent the most likely coordinates of the fall-of-shot. 


FIGURE 3. TYPICAL CORRELATION FUNCTION 
FOR CONSTANT FIRING RATE 


3. COMPUTING THE MISS DISTANCE STATISTICS FROM FIRING DATA 


Volume II of this report proposes a program to derive the statistical 
error parameters primarily trom the theoretical considerations of system 
performance. It has been the practice to measure or compute these values 
from the results of firing trials. Since a computation based on a finite 
amount of such firing data yields only an estimate of the quantity under 
consideration, one must distinguish between this estimate and the true value 
of the parameter. We will identify the former by placing a circumflex over 
the appropriate symbol. Thus, for example, i is the best estimate, on the 


basis of firing data, of the true average aim point u. 


The objective of this section is to indicate the procedure to be followed 
in arriving at estimates of s, u, a,f, and d from records of miss-distance 
data. 


It is necessary to differentiate between the case where a single record 
is available and where many engagements between the same target and gun are 
recorded. In the first instance, we can only compute delivery accuracy for a 
Single occasion, whereas in the latter, we may in addition obtain averages over 


many cngagements. 


Assume then that tnere is avarlable the miss-distance record of a single 
engagement during which n shots were fired. Figure 4 shows a typical example. 
The first step consists of computing an estimate of the true average aim point, 


wnich is given by 


ae 
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Figure § MISS DISTANCE RECORD 


Uver the long run, these estimates will yield the smallest error and also 
have the highest probability of being equal to the true aim point. In more 


traditional terminology, a is the mean point of impact (MPI). 


In general, the MPI computed on the basis of the n impact points will 
differ from the true MPI by a random quantity called the MPI error. The 


standard deviation of this randomness is given by 


S.D. of the MPI estimate = S.D. of impact points 
von 


Thus, the MPI estimate actually approaches the true value of the MPI as the 
number of impact points used in the computation becomes large. Note that this 


formula is only valid when the shots are fired independently of each other. 


With only a single record available, it is impossible to estimate §. This 
é quantity must, therefore, be obtained from other considerations. If there were 
m> 1 such records, however, the variance of the MPI's could be estimated by 
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where uo andu_— are the MPI's for the q”” record. 
*q %q 
The next quantity to be estimated is the total dispersion in the rounds. 
This consists not only of ballistic dispersion, but includes also the effects 
of the time-varying errors on the fall-of-shot. Under our previous assumptions 


on the crror distributions, the total dispersions in the x and y direction are 
a) 4 + ? 
given by de + ay and q., + a. llere, d accounts for the baliistic dispersion 


and a is the contribution of the time varying crrors. An estimate of total 


ORDER ree 
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dispersion may be obtained from a single firing record with MPI, a, according 
to the formula 


7S ; an , 2 
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Since the ballistic dispersion, d, could be derived for a given operational 
situation by combining the increments contributing to this dispersion in 
accordance with techniques to be developed in the program of Vol. II, the last 
formulas can be solved for a. 


It remains to estimate the correlation functions. ‘f 7 and fy from the 
firing trial data. By definition, 


feos fy om #1 


A good estimate at the other values of the argument J may be obtained from 
the equations 
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Where j is the number of intervals separating the rounds whose correlation is 
being estimated, n is the total number of rounds fired, and (7-/) is the 
number of sample data points separated by exactly j intervals. 


The equivalent constant correlation r referred to in the previous section 


is defined by the formulas 


The above is one of several approximations which could be employed, and is 
cited in Reference 3 as providing reasonable estimates, based on quadratic 
weichting of the Ji) . In this expression, m-1 is the number of different 
coefficients Pu which are considered in arrivinn at the estimate. 


4, DELIVERY ACCURACY 


Weapon effectiveness is a measure of the weapon system's ability to destroy 
a target. As such, its value is not only a function of the gun and associated 
fire control equipment, but to a large extent depends on the type of target 
under attack and the warhead lethality. Thus, weapon effectiveness is a conven- 
tient parameter for the purpose of mission planning, where the targets are well 
defined, but suffers serious disadvantages in the comparison of existing and 
planned systems. Delivery accuracy avoids some of these shortcomings by measuring 
only the system's ability to place the warheads within a specified distance of 
the target, regardless of the nature of the target and warhead. 


Because of the random variations of time-varying and round-to-round disper- 
sion errors, the impact point of any particular round is, for a given engagement, 
random. This implies that delivery accuracy must be measured in terms of the 
statistical quantities which specify the distribution of impact points. Ona 
single engagement these are the previously defined parameters u, a, and d. 

To this list must be added 5S, the dispersion of the systematic error when the 
delivery accuracy over many occasions is of interest. 


It is possible to accept these parameters as defining delivery accuracy. 
llowever, in most situations the quantities of major interest are not the miss 
distance statistics but the likelihood of achieving a certain number of hits, k, 
within a well-defined area about the target. Thus, we are led to make the 
following , 


Definition: Delivery accuracy is the probability distribution of the nurber of 
hits on a fixed area containing the target, i.e., the probability of exactly k 
hits when n rounds are fired, for all O¢ k Cn. 


If now the target is of such a nature that the damage inflicted on it is 
only dependent on the number of hits obtained on an area about it, the actual 
damage during an engagement may be assessed probabilistically by taking the product 
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of the "damage function" and the hit distribution. For example, if Q(k) is the 
probability that the target is killed when k hits are received, the unconditional 
probability of a kill in an engagement during which n rounds are fired is 


hap 5 CA) Ly) 


k?0 


where PCa) is the probability of k hits in n rounds. 


In situations where the damage depends on the actual coordinates of the 
impact point, a multiple-area equivalent target may be defined, to which similar 
considerations as above then apply. 


Many targets have the property that they can be incapacitated with a single 
hit in the vulnerable portions. If the function Q({k) refers to just that part of 
the total target area then 


O;% =O 
Pia) = 
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Substituting into the previous forma we have then 


nw 
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which shows that the target kill probability is equal to the probability of at 
least one hit, i.e., the cumulative hit probability. It is only in such cases 
as described here that the istter parameter is sufficient for a complete assess- 
ment of weapon effectiveness. Fortunately, these include many situations of 
interest. 


Although the hit distribution is derived from the statistics of the miss 
distances, the relationship is by no means obvious. In the remaining part of 
the present section we elaborate on this point. 
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Define Hy to be the event that the i round results in a hit, and P(Hi, Hi,.+. 
Hi,) as the probability that rounds iy iy, eas i, are hits. Then the probability 
of exactly k hits when n rounds are fired is given by 


2e-k , : 
2) S00) by, 
2=0 


where, by definition, 


Ey= 2 Ply hy Mi Zt 


by Seg <b 


The summation in the last equation is over all Q) sets of k integers chosen 
from the set (1, 2, ... n). 


The cumulative distribution function, i.e., the probability of at least 
k hits out of n rounds, is denoted by R, tk) and given in terms of the E. by 


wh, ry hk-Lsi 
R, lh) =) (tA JE 

According to the above fais the hit distribution is completely determined 
by the probabilities of hitting with a given subsequence of rounds, for all such 
subsequences. Our next task is to relate the miss-distance statistics to these 
quantities. The actual computation involves certain multi-dimensional integrals 
over finite limits, a problem for which no closed-form solution exists. An 
approximate solution, however, may easily be obtained by replacing the correlation 
function © by its constant equivalent r, providing a weighting function for the 
integrand and extending the limits of integration to infinity. The details of 
this method may be found in Reference 1. Here we restrict ourselves to stating 
the result which is given by 
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In these equations C. are Cc, are parameters depending on target shape and size. 


For a circular target of area A, for example, C. = cy “/Az. For a rectangular 


target of sides 2w and 2h, C. = wilt .C = nie For a general area A, 

the relation C. C = Al7 holds. The matrices 1 ‘ and A y are k by k 
iden cays : 2 2 2 2 2 2 

whose main-diagonal entries are d. +s + 8 and d, + sy + a, , and whose 


; : 2 2 2 2 , 
off-diagonal entries are Ss. +a fT, and sy + ay Ty: respectively. The 


Symbol I denotes the k by k unit matrix and Uy and U, are k x 1 column matrices 


whose entries all equal u and Uy» respectively. 


The matrices in the above formulas are simple enough so that the determinants 
and inverses may be evaluated in a straight-forward manner, The result, for 


either the x or y component, is 


[a rri/=plerraldrat raiar(srniyferdrarre 
and 
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Note that, as a result of replacing the correlation function by the equivalent 
constant value r, the above terms are independent of the round indices iy, io Size 
i, Consequently, the probability of hitting with a certain subsequence of rounds 


is the same for all sequences of k rounds. 


It now becomes necessary to distinguish between two cases: the hit distribution 
on a single engagment for which the systematic error is u; and the average hit 
distribution over many such engagements, where the average systematic error is 
zero but varies from one engagement to the next with standard deviation, 5,In the 
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former we set S = 0 in our general formula for hit distribution, whereas in 
the latter u = 0. Thus, we obtain the following equations: 


A. Hit distribution on a single engagment with systematic error u, time-varying 
error standard deviation a, equivalent correlation r, and round-to-round 
dispersion error standard deviation d. 
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B. Hit distribution over many engagements with systematic error standard 
deviations,s, time-varying error standard deviation, a,and equivalent correlation, 
r, and round-to-round dispersion error standard deviation,d. 
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Two special cases merit mentioning. If the rate of fire is fast enough 
so that the aimpoint changes little during the engagement, we have approximately 


ey eee. for all 7 


and consequently, rs r = 1. This firing mode is variously described as 
burst or shotgun fire. 


If on the other hand, the rate of variation of the time-varying errors 
exceeds the firing rate, the correlation in these errors becomes zero and 


r* Ty = 0. The impact points of the rounds become statistically independent, 


and the hit distribution on one engagement reduces to 
BCA) = (BYP tu) J [1- Poy)” 


where P(H) is the single-shot hit probability, given by 


Cy Cy 


PI) * Loe +2(dhirae]] [Cy +2 (dy 4 ay*)/ 


2 2 
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5. SINGLE SHOT HIT PROBABILITY, CUMULATIVE HIT PROBABILITY, 
AND EXPECTED NUMBER OF HITS 


Some of the most useful quantities derived from delivery accuracy, and at 
the same time, the simplest from a computational point of view, are single shot 
hit probability, the probability of at least one hit in n rounds, and the expect- 
ed number of hits. If the target can he destroyed with one hit, then the former 
two are the only measures which need to be considered. In the following, we 


present the formulas as they specialize to these three cases. 


I, Single Shot Hit Probability 


The probability of hitting the target with a single shot, regardless 
of the results of ali other shots, is obtained by substituting into the 
equations in the previous section the values n = k = 1. For case A (single 
engagement) we obtain, after some manipulation, the result 


Cx Cy 


P(H) = Vie r2(aer ar) {Cj til y ray )/ 


2 
Ay 


ee Oy 42a," + az ) = G oaiiee przaia)} 


When the target is square or circular of area A and the error statistics 


are cqual in both coordinates, this reduces to the simle form 
27 ue 


A 7 
Sar ay SOLE e  Atandya™) 


which is equivalent to the well-known Carlton formula. 
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For case B (multi-engagements) the single shot hit probability becomes 


F(A) = OE RR SR 
LGf te (Ext + apres] [Cf 42 (Sy + ay +ay*)] 


Again, for square or circular targets of area A a-d equal error statistics, 


this simplifies to 


Pe eae. eaten 
PU) A + 20(5*4a*1a") 


II. Cumulative Hit Probability 


The probability of hitting with at least one round (k=l) when n rounds 
are fired is given by 


4 e 
ent 
LK ,, (4) = > (-°~ £, 
n é 
c=Z 
where E. is the probability of hitting with i rounds.* For case A we, there- 


fore have 
Ritz og ey cA a 
fet [tsp dla dy ty) +213 a8] [ig eae vay Tay] 
Cy! 
x Vl Gj +2laj ry yagt)e2iry 43] /CFreldy dty' 7, Ay) PL 
aye 


Ee ITT DIT PTT, 


X CAP y, = 
is 2 2 a 2 : z 7 
LS + £(ay tay" Fy ay )4264, a3] 
“Observe that the combinatorial tera (Klee) which appeared in the original 
formulation for RA (k) is unity for K = 1. 
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When the target is square or circular of area A and the error statistics 


are equal in both coordinates, this reduces to the simpler form 


A 


ft 
GS) gl (5 ee A eS A cern ec Aine en 
Fn - Cu (2 At Te 4a*-Ya*) - / ae ‘a*-ra*) +27 ra* 


| aa LH2« 
xX €Xp- FA A+2M(d +a*-ra*)+2Nira* yh 


For case B the corresponding formulas are 


a. 0. Z 
LM 
fy" lee ztaped}-nag)+2ilgf+ 5 4 [of aldirad aby 


Cy +2 (ay ray . yay y wae (Syrray icy" Yay 14/5; 4;') ef 
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Ky * Jf “L) (0 licarisam +774) vA os A iF rain / 
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III. Expected Number of Hits 


This measure is simply given by the product of the numper of rounds fired 


and the probability of hitting with a single shot. Thus 


El(H) =n2P{H) 


where ) is the single shot hit probability appropriate to the case under 


consideration. 


It should be observed that this expression is valid only under circum- 
stances where the single shot hit prohability does not change from round-to- 
round. This requirement does not necessarily imly that successive rounds must 


be independent, (uncorrelated) however, as illustrated by a simple example. 


Consider a burst of two (2) rounds with the following hit probabilities. 
Denote the event "a hit on the jth round" as He and the event "a miss on the 
i? round" as H;. 

Round 1: P (Hit on first round) = P (Hy) = 0.60 
P (Miss on first round) = P qi) = 0.40 


Round 2: Conditional Probabilities 
P (Hit on 2nd round | Hit on first) = P (H, | I) = 0.7 
P (Miss on 2nd round| Hit on first) » P (HZ | 1.) = 0.3 
P (Hit on 2nd round | Miss on first) = P (li, | ii}) = 0.45 
P (Miss on 2nd round | Miss on first) = P (ii, | Th) = 0.55 


The first and second rounds in the above example are correlated, since 
the probability of a hit on the 2nd round depends on whether or not a hit 
occurred on the first round. However, observe that the single shot hit 
probability is the same for both rounds. That is, 


P (iy) = 0.60 
and 


P(N) =P (Hy Hy) + PH) + P| AD PD 


= (0.7) x (0.6) + (0.45) (0.4) 
= 0.60 


2) 


-25- 


In this case, the expected number of hits E (H) is n- P (H) = 2x .6 = 1.2. 


In the more general case, where the rounds are correlated but the single 
shot hit probability does not remain constant over the salvo, the expected 
number of hits may be obtained by subdividing the salvo into intervals in which 
the single shot hit probability is essentially constant, and summing the number 
of expected hits over each interval. For example, consider a very simple case 
where n = 2 rounds and the first and second round conditional probabilities 
are as follows: 


P (i,) = 0.6 
P (H)) = 0.4 
P My | H)) = 0.8 
P (i, | Hy) = 0.2 
P (H, | Hy) = 0.6 
P (ii, | Hj) = 0.4 


Again, the first and second rounds are clearly correlated; but now, 
P (Hy) = 0.6, 
while 
SS 2 
P (Ht) 0.72. 
In this case, the expected number of hits is given by 
E (tl) = n° P (ty) Nn, P (H,) 


Where ny and ny are the number of rounds in the first and second intervals, 


respectively, over which the probabilities are constant (llere, ny "mn, * 1). 


Thus, E (i) = 0.6 ¢ 0.72 = 1.32. 
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6. A DETAILED EXAMPLE 
When employing this manual, it is expected that relationships permitting 
military planners to derive values of the parameters, a, d, s, u, and J for 
specific situations would be available from other sources. The values of X , 
£2, and ff necessary to enter the appended tables could then be computed. The 
program of Volume II is intended to provide these relationships. 


The example cited here, in addition to computing hit probability, demon- 
strates how values of the statistical parameters relevant to the firing of a 
specific salvo might also be derived from the results of the firing; provided an 
applicable value of round-to-round dispersion is available. An examination such 
as that presented here might be performed on the results of an actual firing 
exercise in order to verify that pretest estimates of the parameters are correct 
and to provide a basis for preparing improved estimates of these parameters and 
hit probabilities for future logistic and fire planning. 


Let us consider the case of a gun firing twenty rounds at a square target 
30-meters on a side and compute the hit probabilities for both single and multiple 
engagements. Let us also assume that a previous firing exercise resulted in the 
pattern of impact points shown in Figure S and that from other considerations, 
it is known that the standard deviations of the round-to-round dispersion and 
systematic errors are 10 meters and 30 meters respectively, equal in both coordi- 
nates. The impact point data is presented in Table 2. We assume that all spotting 
corrections have been subtracted from the impact points and that the resulting 
sequence of shots is statistically stationary. 


Round Number (k) 


EERLEEEBECEEEEEEEoE 


X-Coordinate of h , Impact Point (X,) 


Y-Coordinate of k*" Impact Point (Y,) 


TABLE 2, FIRING DATA 
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Figure S SAMPLE FIRING TEST RESULTS 
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The first objective is to compute an estimate of the true average aim point u 
for the engagement resulting in the above impact point data. Frem the equations 


oe Z rae Pa ay 

Pe ae 3 P oS ag ee 

7 a 2 ae, ge Re 2 4 
rae d Z=Z 


we obtain 


fy "SEMETERS 3 hy = EF METERS 


Thus, the location of the mean point of impact is 58 meters to the right =-nd 
65 meters beyond the target center, as indicated in Figure 5S. 


The next .tarcity to be estimated is the standard deviation of the time- 


varying error. Since the round-to-round dispersion d is known to equal 10 
meters in both coordinates, we can use the firing data of Table 2 and obtain the 


total dispersion, 


| [ole + a2 = Vax (X; - fey)” +28 METERS 
cf 
a = 2 
dy + ay = ES" (Y- 4y)” = 2B METERS 
Les 


and from that, the time-varying error dispersion, 


Ay ~2LMETERS 5 By * 26 METERS 
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Thus, if the firing which resulted in the data of Figure 5 had been continued 
indefinitely, approximately 40 percent of the rounds would have fallen in a 
circle of radius 28 meters, centered at the mean point of impact. In the actual 


firing 50 percent of the rounds were so located, as shown in Figure 5S. 


The correlation between rounds 7 units apart in the sequence can be estimated 


from the formula 


Bee ee Z 4 
at ali) = Eid” [> ae] [Kyo 4) 
t=L 


and the corresponding formula for the y-component. In this expression, j is the 
number of intervals separating the rounds whose correlation is being computed, 
n is the totai number of rounds fired, and (n - j) is the number of sample points 


for vshich data are available. 


For example, we have from Table 2, and the above values for a, cd, and u, 


with j = 10 and n = 20 rounds: 


5 r 
Px et OMeist Baga [Oy nd hn hn) + Oa Ue Vg Ud 
# UX, - Ux) X_- Uy) -. 2. 


+ (X,. — Uy MXQ~- uJ 


p (10)* /Zpy2 jé i (2)\-2¢) + (-4N-48) + (-L2Y- 68) + (12)f-39) 
* (20-26) +(22N30) +(f2N32) + 2ENG2Z) 
¢(20NI0) + (-8-e)] /} 


= 0.39 
The first ten values of Sf. and f, ave tabulated in Table 3. Note that 


the correlation between rounds separated by 10 intervening rounds is still 


not negligible. 


Correlation In 
X-Corrdinate 


P,9) 


Correlation In 
Y-Coordinate 


Ay 


TABLE 3. FIRING DATA CORRELATION 
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The accuracy of the correlation estimates is reasonably high, considering 
the small amount of data available, for small values of aa , but deteriorates 
rapidly as J increases beyond 10. This is, of course, due to the fact that 
the number of terms which are averaged decreases linearly with J . For example, 
the estimate of FP (19) is based on a single term, whereas the formula for 
Pa) is an average of 19 terms. This is the reason why the data in Table 3 
extends only to / = 10. 


With the aid of the formula 


yre-7 


2: BID) 
ey 5 einen 


nr 
a (mJ)? 
FL 


and J#%= 10, the equivalent correlation r becomes, approximately, 
Py =hy = OLZ 
Sees | 4 
Here 77-Jis the number of terms from which a good estimate of P can be made. 


This concludes the computation of the miss-distance statistics. 


Even though the magnitude of r is quite small in this example case, a 
complete set of computations is now outlined, including the correlation parameters, 
in order to demonstrate the methodology presented in this manual. 
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Values for the probability of at least one hit on a target when n rounds 
are fired, RL (1), and for the single shot hit probability, P (H) = Ry (1), 
may be obtained from the tables appended to this manual. These tables have 


been prepared in generalized form in terms of the following composite 


parameters. 


I, Single Engagement Case: 
2 
Eee d+ a*(1-y) 
A 


se PE 
p A 

2 
aa 


II. Multiple Engagement Case: 


Oh Kee d+2(1-F7) 


=o 


Standard deviation of systematic error, many engagements 


z 
z 
e 
te J 
e 


Actual single engagement systematic error 
Standard deviation, round-to-round dispersion 
Standard deviation, time varying error 
Correlation function for time varying error, ca 
Equivalent average value of 
Number of rounds 
Number of hits per salvo 

- Area of a nearly-square or circular target 


>eX¥ YHSAARG 
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Since the target in the above example is square, and the miss distance 
statistics are approximately equal in the two coordinates, we may obtain the 
desired probabilities from these tables. The parameter values developed in 


the above example are as follows: 


S = 30 meters (Standard deviation of systematic error; assumed) 

u = 62 meters (Actual systematic error; computed from data of 
Table 2) 

d = 10 meters (Dispersion error; assumed known from previous 
firing) 

a = 26 meters (Time varying error; computed from knowledge of d 
and data of Table 2) 

r= .12 (Equivalent corelation parameter; estimated above) 

n= 20 (Number of rounds fired in salvo) 

k=] (Number of hits per salvo) 

A = 900 meter” (Assumed 30 x 30 meter target) 


Using these values, we compute: 


. Pe a eee?) 


a = 43 (Single Engagement) 
O 


x a OF 

A 

2 
ve _ as = 2. OF (Single Engagement) 

2 2 
=o ~- 1.09 (Multiple Engagement) 
, 2 

| 
| | 


(Multiple Engagement) 


262 


Soaeanere geese 


ci ce en me 
ec ance AAR AOE a OA CIID AOS SAE AEB . 


eee 


After some interpolation in the appropriate tables or graphs, we finally 
obtain for single shot hit probability, cumulative hit probability, and 


expected number of hits the values 
P(H) = 0.0025 
R,, (1) = 0.0477 
FE (H) = 0.050 


over the single 20-round engagement. 


If for purposes of this example, 


one makes the unrealistic assumption 
that d, a, and r remain constant over 


a number of engagements, 


P(H) = 0.0795 
Roo (1) = 0.687 
E(H) = 159 


as averages over many engagements. Table 4 summarizes the example results. 


In all of the above examples, the expected number of hits E (H) has been 
computed from the expression E (H) =n ¢ P (Il), under the assumption that 
P (H) does not vary significantly from round-to-round. 
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If one chose to compute the probabilities P (H) and Ro (1) ignoring the 
correlation effect, the following procedure could be employed. The total 
2 
d' +a 


dispersion, d', was defined previously as d' = 


2 = 28 meters, where 


d and a are the round-to-round dispersion and the time-varying error dispersion 


components, respectively. If the correlation is to be ignored, one can set 


r= a= 0, and, in effect, ascribe the total observed dispersion d' to the 


round-to-round dispersion component. 


Now, for the single-engagement situation, the parameters have the values: 


Bp -: 
(: 


0.87 


0 (0.01 was used in entering the appended 
tables) 


4.3 


and the resulting values of P (lH), Rog) and E (H) are: 


P (H) = 0.0026 

Rag (2) = 0.0518 

E (H) z 0.052 

° Similarly, for the multiple engagem-nt case, assuming S = 30 meters, we 
have: 

Ow = 0.87 
poo 
| 


and the resulting v-lues of P (il), R55 (4), and E (H) are: 


P(li) 
Rog (L)® 
EQ) 


0.0786 
0.707 
1.57 


-38- 


Observe that these values do not differ markedly from those computed 
under the previous set of assumptions (namely, with a = 26, r = 0,12, and 
d = 10), as would be expected with such a small value for the correlation 


found in this particular example. 
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7. TABLES §& GRAPHS 


The tables appended to this report provide computed values of the prob- 
ability of at least one hit on a target when n rounds are fired, RQ, (1), for 
nearly square or circular targets of area A and equal miss distance statistics 
in both coordinates. They include as a special case the single shot hit 
probability Ri (1) = P (i). 


The tables are used by computing values for the generalized parameters 
x, B & a which are defined in the previous section as functions of the 
target and error parameters A, a, d, r, u, and S. Intermediate values may 
be obtained by interpolation as illustrated in the above example. Values of 
RA (1) are tabulated for any number of rounds n between 1 and 30 and values 


of the remaining parameters are as follows: 


aH = 0.01, 0.20 (0.20) 1.00 ©.50) 10.00 (1.00) 16.00 
J = 0.01, 0.20 (0.20) 1.00 (0.50) 10.00 (1.00) 16.00 
y¥ = 01) S, 10 


Figures 6 through 17 are plots of cumilative hit probability, RA (1), 
for 10, 20, and 30 rounds fired and selected values of Y. They are included 
primarily to demonstrate the trends resulting from changes in the several 
input functions. A constant ordinate or vertical scale has been used throughout 
in order to permit direct visual comparison of hit probability. The abscissa 
scale, however, has been increased for those curves based on larger values of 
Y in order to clarify the overlapping situation which exists at the lower 


values of GX. 
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Appendix 
TABLE OF CUMULATIVE HIT PROBABILITY 
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